Peak effect and dynamic melting of vortex matter in NbSe 2 crystals 
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We present a mode locking (ML) phenomenon of vortex matter observed around the peak effect 
regime of 2H-NbSe2 pure single crystals. The ML features allow us not only to trace how the shear 
rigidity of driven vortices persists on approaching the second critical field, but also to demonstrate a 
dynamic melting transition of driven vortices at a given velocity. We observe the velocity dependent 
melting signatures in the peak effect regime, which reveal a crossover between the disorder-induced 
transition at small velocity and the thermally induced transition at large velocity. This uncovers 
the relationship between the peak effect and the thermal melting. 

PACS numbers: 74.25.Qt, 74.25.Sv 
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Discontinuous jumps in equilibrium magnetization ob- 
served well below the mean field line in clean high T c 
cuprate superconductors [1] , have been widely recognized 
as a hallmark of the thermodynamic melting transition 
(MT), which separates a vortex solid state, where elas- 
tic interaction dominates and quasi-long ranged, crys- 
talline correlations develop in vortex structure, from a 
vortex liquid state in which thermal fluctuations disrupt 
the crystalline order and the shear rigidity vanishes. 

The melting signature is often accompanied by the pro- 
nounced peak anomaly of the magnetization or critical 
current [2], known as the peak effect (PE), originating 
from the rapid softening of the vortex lattice and the 
random pinning potential due to disorder quenched in 
a host material. The close proximity of the MT on the 
PE has led to a reexamination of the physical properties 
close to the second critical field H c2 , especially on low T c 
materials like Nb [3-5] and NbSe2 [6,7]. Various exper- 
imental results evidence the pinning-induced structural 
transformation into a disordered array of vortices, indi- 
cating the dominant influence of the quenched disorder in 
the PE regime [3,7-9]. However, the presence of the MT 
and its relation with the PE remain controversial [4,5]. 

The effect of the quenched disorder changes remark- 
ably when vortices are driven by a transport current. 
As proposed theoretically [10-12] and experimentally 
[6,13,14], at low drive the "shaking action" of the random 
pinning due to the quenched disorder disrupts largely the 
internal periodicity in moving structure of driven vor- 
tices, while at high drive the influence of the disorder di- 
minishes and the elastic interaction becomes dominant, 
resulting in moving solid states at large velocity. 

This leads to an unique opportunity to study whether 
the "driven lattice" undergoes thermally induced MT and 
also whether the transition occurs in the disorder domi- 
nant, PE regime. In this letter, we present experimental 
evidences for the dynamic melting transition(DMT) of 
driven lattice observed just above the PE of NbSe2 by 
means of mode locking (ML) experiment. ML is a dy- 
namic synchronization between rf drive superimposed on 
dc drive and collective lattice (elastic) modes excited over 
driven vortices at internal frequencies fi n t = qv/a with 



the average velocity v, the lattice spacing a and an integer 
q [15,16]. When the internal frequency and the external 
drive frequency are harmonically related, i.e. fi nt = pf 
or v = (p/q)fa with an integer p, the elastic modes are 
dynamically locked on rf drive. We observe this dynamic 
resonance as multiple current steps in a current-voltage 
(IV) curve or multiple peaks in a differential conduc- 
tance (dl/dV) curve. The ML features allow us not only 
to trace how the shear rigidity of driven vortices persists 
on approaching H c2 , but also to demonstrate the DMT 
of driven vortices at a given velocity. We will show that 
the velocity dependence of dynamic melting signatures 
unveils the relationship between the PE and the MT. 

Platelets of 2H-NbSe2 pure single crystal used in this 
study were grown by an iodine vapor transport method 
[17]. The thinnest three crystals (typically 10/im thick) 
were cut in the bar shape and cleaved with no signif- 
icant optical surface damage. Contacts for the four- 
probe method were made by indium solder. The crystals 
have the superconducting transition temperature T c of 
7.2 K determined from the mid-point of resistive tran- 
sition. The transition width is about 50 mK between 
10% and 90% of the normal state resistance just above 
T c (e.g., R(7.4K)). Typical values of residual resistance 
ratio i?(295K)/i?(7.4K) for the crystals are ranged be- 
tween 30 and 40. However, no significant change of the 
room temperature resistivity po of 1.2 x 10 -6 fhri was ob- 
served in crystals from different batches [17]. This allows 
us to estimate the thickness of the cleaved crystals from 
the room temperature resistance i?(295K). 

The measurements discussed in this study were carried 
out on a cleaved crystal with a dimension of 0.72mm(Z)x 
0.63mm(w)x 0.9(i)/im. Figure 1(b) shows results of H c2 
on the crystal obtained from a measurement of flux flow 
resistance with magnetic field applied perpendicularly to 
the platelet. As displayed in the inset to Fig. 1(a), H c2 
is determined from the intersectional point between the 
linear extrapolation of flow resistance indicated by a solid 
line and the normal resistance by a dotted line [18]. The 
slope of H c2 on T at T c (= -S) is -0.72 T/K, in good 
agreement with magnetization results [7,17,19]. Shown 
in Fig. 1(a) is a plot of the critical current I c vs mag- 
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netic field H, obtained by ramping field up (and also 
down) after cooling the crystal to 4.2 K in the zero field 
(indicated by arrows in Fig. 1(b)) [20]. Here, I c is deter- 
mined from IV curves by a 10fj,V/m criterion. It exhibits 
a peak behavior below fioH c2 of 2.14 T, characterized by 
the peak field H p of 1.93 T. We note that, because of 
inhomogeneous edges in the strip geometry, I c below H p 
has both bulk and edge contributions [21,22], resulting 
in broadening for the lower part of the PE. 

Figure 2(a) shows a series of dl/dV — V curves mea- 
sured in 1.90 T (just below H p ) with superimposing 3 
MHz rf current of various amplitudes [23]. Applica- 
tion of a relatively large rf current of 16.7 mA, for in- 
stance, induces clear conductance ML peaks at equidis- 
tant voltages. First peak corresponds to the fundamental 
{p/q = 1/1) an( i others are higher harmonics (2/1, and 
3/1). Small peaks at sub- harmonic conditions with q > 1 
are also detected. As expected in the velocity condition 
(v = (p/q)fa), those resonance voltages (oc v) increase 
linearly with rf drive frequency (see Fig. 2(b)). Thus, 
by changing /, we can detect the shear rigidity of driven 
vortices at various velocities. 

In Fig. 2(c), we show how the fundamental resonance 
depends on rf current I r j. Here, the fundamental ML 
current width AI is obtained by integrating the conduc- 
tance ML peak with respect to the flux-flow-base line 
(see the inset to Fig. 2(c)). It shows an oscillatory be- 
havior with I r f , which agrees qualitatively with squared 
Bessel function of the first kind displayed by a broken 
curve. This behavior is expected when the random pin- 
ning due to the quenched disorder excites elastic modes 
in driven clastic lattices [15]. In other frequencies and 
fields, similar oscillatory behavior of AI on I r f are ob- 
served in different rf current scales. For simplicity, we 
focus mainly on the first maximum AI max of the funda- 
mental ML width in the following analysis. 

Next we turn to the magnetic field dependence of the 
ML resonance. In Fig. 3(a), we show how the fundamen- 
tal ML voltage Vi/i evolves with H. As observed, Vi/i 
for 7 MHz increases with H with downward curvature. 
This is in good quantitative agreement with the funda- 
mental ML-voltage condition displayed by a solid curve; 
Vi/i = Af&ol/ao with the equilibrium lattice spacing 

a (= y/AfyjB), a factor A(= 2/y/S), the sample length 
l(=0.72 mm), and the vortex density B = hqH. Note 
that no fitting parameter is used for this comparison. 

As magnetic field approaches /Uo-ff C 2(=2.14 T), how- 
ever, the ML resonance suddenly disappears at a certain 
field. An example is displayed in Fig. 3(b), where the 
dl/dV curves measured with superimposing a 7 MHz rf 
current in different fields are shown. The fundamental 
ML peak is clearly observed in lower fields. As H is 
ramped up, the peak becomes small and it seems to van- 
ish around /x o i?=2.01 T, above which the dl/dV curves 
become featureless. 

The characteristics are better displayed by plotting 
AI m ax against field, which is given in Fig. 4(a). AI max 



for 3 MHz, for instance, starts to show a rapid drop 
around 1.98T (defined as the onset field H on ) and dis- 
appears around 2.00T, while it is insensitive to H below 
Hon- We define the vanishing field H* from the linear 
extrapolation to AI max = as indicated by a solid line. 
Above this field, no ML resonance appears at any ampli- 
tude of I r f , indicating the absence of the shear rigidity 
in driven vortices. Thus, the vanishing field should mark 
the DMT from coherent to liquid like incoherent flow 
states [15,16]. We note that this phenomenon emerges 
just above the PE where I c drops, and therefore it should 
be driven by the thermal fluctuations, in addition to the 
quenched disorder. 

The reduction of the disorder effect is expected in the 
ML features measured at high frequencies (large veloci- 
ties). As shown in Fig. 4(a), AI max measured at a high 
frequency of 9 MHz exhibits a sharper drop at slightly 
higher fields. Namely, both the onset and the vanishing 
fields increase slightly, while the field range where the 
linear drop of the ML width occurs becomes narrow. We 
observe even sharper phenomenon at high frequencies, 
indicating the sudden disappearance of the shear rigidity 
in driven vortices at large velocities. 

The striking frequency (velocity) dependence is also 
observed in the onset and vanishing fields, characteriz- 
ing the dynamic melting signatures. In Fig. 4(b) results 
of H* and H on are plotted against velocity determined 
from the fundamental velocity condition of v = fao. As 
observed, H* (open symbols) increases with v and ex- 
hibits a saturation behavior toward w2.01 T at larger 
velocities. H on (solid symbols) also exhibits a similar 
saturation behavior toward the nearly identical field at 
larger velocities. The difference between H* and H on at 
a constant velocity decreases with increasing v and be- 
comes negligible experimentally at large velocity. These 
behaviors imply that at small velocity the influence of 
the disorder facilitates the occurrence of the DMT, while 
at large velocity it becomes negligible and the thermal 
effect dominates on the phenomenon. Thus, the velocity 
dependence of the dynamic melting signatures reveals a 
crossover between the disorder- and the thermally-driven 
transitions. 

Regarding the vanishing field at a given velocity as 
a crystallization velocity v c at the field, the saturation 
behavior is qualitatively similar to the dynamic order- 
ing picture proposed by Koshelev and Vinokur [10], in 
which combined influences of the thermal fluctuations 
and the quenched disorder on vortex dynamics near the 
thermodynamic melting transition arc taken into ac- 
count. In this model, the crystallization velocity di- 
verges as v c = vo/(l — T/Tm) on approaching the melt- 
ing point Tm from below. As displayed by a solid curve, 
the H* data are well approximated by a similar func- 
tion of u c = vo/(l — H/Hm)" with the melting field 
Mo-Hm = 2.01T (indicated by a broken line), a ~1 and 
t>o = 0.33mm/s over nearly two decades of velocity [24]. 
This allow us to identify the saturation field as the ther- 
modynamic melting point for the vortex lattice. Thus, 
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the sharp drop of ML by the saturation field signifies the 
thermally induced, DMT of driven lattice. 

Let us compare the saturation field to a recent quan- 
titative theoretical result for the thermodynamic melt- 
ing line proposed by Li and Rosenstein [25]. The MT is 
set in whe n the l owest Landau-level-scaled temperature 
cit = -(y / N Gl /2Trth)- 2 / 3 (l -t-h) = -9.5, with their 
Ginsburg number N Gi = {{k B T cl / '4tt 2 // ff c (0) 2 f 3 ) 2 ) /2, 
a reduced temperature t = T/T c , and a reduced field 
h = H/H c2 (0). Using the anisotropy parameter 7 = 3 
[19], the Ginzburg-Landau (GL) thermodynamic field at 
T =0 n H c (0) = 0.24 T [19], and the GL second critical 
field n H c2 (0) = ST C = 5.18 T for the present crystal, 
we find that the saturation field of 2.01 T at 4.2 K gives 
ar =-10. A nearly similar result of a,T=-9 is also found 
in other crystal with /zo-Hm=2.02 T at 4.2 K. These are 
in good agreement with the theoretical result. 

Next, we discuss flow states near the DMT. In the 
regime between the onset and vanishing fields, an anoma- 
lous dependence of AI on I r f is observed, of which exam- 
ple is shown in Fig. 2(d). There is a minimum rf current 
to observe ML, indicating the influence of the quenched 
disorder disrupting the shear rigidity in driven vortices. 
Such threshold behavior in ML favors with a smectic flow 
(or a layered liquid) state [26,27], in which only the peri- 
odicity perpendicular to the flow direction is preserved. 
Thus, the H* curve would mark a dynamic transition be- 
tween the smectic and the liquid flow states. Meanwhile, 
for H < H on , the threshold rf current for ML disappears 
and AI shows the the Bessel function like oscillatory be- 
havior on I r f, suggesting a moving solid state. Thus, the 
solid-smectic transition would occur at the H on curve. 
Similar two step process of the DMT has been proposed 
theoretically [11,12]. 

Finally, we compare the melting curve of H* to the 
peak behavior of I c . As indicated by a dotted line in 
Fig. 1(a), the dynamic melting point (or Hm) at large 
velocity emerges in the higher part of the PE. On reduc- 
ing velocity, H* approaches H p and seems to disappear 
(or be below our experimental resolutions) around 1.9 T 
just below H p . A recent study on edge and bulk trans- 
port [22] reveals that the bulk component of I c begins to 
show a peak anomaly just below the peak temperature 
of I c . This prompts us to relate the disappearance of 
H* to the disordering transition for the PE. Thus, the 
H * curve would be a locus of the melting point connect- 
ing smoothly the PE with the MT as a function of the 
disorder effect (velocity). 

In summary, employing ML techniques, we have pre- 
sented experimental evidences for the DMT of driven vor- 
tex lattice observed just above the PE of 2H-NbSe2 single 
crystals. At small velocity the transition driven predom- 
inantly by disorder occurs near H p , while at large ve- 
locity the transition by thermal fluctuations is observed 
by Hm(> Hp). The velocity dependence of the melting 
signatures reveals the smooth crossover between the PE 



and the MT. We hope that this work stimulate 
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FIGURE CAPTION 



Fig. 1 (a) Peak anomaly of critical current I c measured 
by varying magnetic field H at 4.2K. The peak field H p is 
denoted by an arrow. A dotted line marks the dynamic 
melting field measured at high velocity. The inset dis- 
plays how differential resistance Rd measured at a large 
dc current of 30 mA (^> J c ) depends on field. The defi- 
nition of the second critical field H C 2 is given in text, (b) 
Phase diagram of NbSc 2 . Both H c2 and H p are plotted. 
See text for arrows. 

Fig. 2 Mode locking (ML) phenomenon of driven vor- 
tices, (a) Differential conductance (dl/dV)- voltage (V) 
curves measured in 1.90 T with superimposing 3 MHz rf 
current of various amplitudes. For clarity, the origin for 
each curve is shifted vertically. ML conditions are de- 
noted, (b) Frequency / dependence of ML voltages for 
the fundamental and higher harmonics, (c) Rf current 
I r f vs. the fundamental ML current width AJ obtained 
by integrating the fundamental ML peak shown in (a) 
with respect to the flux- flow base line (see the inset). 
See text for a broken curve, (d) AI for 3 MHz vs I r f 
observed in 1.995 T between H* and H on . 

Fig. 3 (a) Magnetic field dependence of the funda- 
mental ML voltage V\j\ measured with superimposing a 
7 MHz rf current. A solid curve represents the funda- 
mental ML- voltage condition (see text), (b) Differential 
conductance curves measured with a 7 MHz rf current in 
fields from 1.96 T to 2.02 T by a constant field step of 
0.01 T. For clarity the origin for each curve is vertically 
shifted. 

Fig. 4 (a) Magnetic field vs. the maximized funda- 
mental ML current width AI max for 3 MHz and 9 MHz. 
For clarity, the origin for 9 MHz data is shifted vertically. 
Onset H on and vanishing H* fields are indicated by ar- 
rows. Their velocity dependencies are given in (b). A 
broken line represents the saturation field. See text for 
solid and dotted curves. 
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